Regular izat ion parameters for the self-force of a scalar particle in 
a general orbit about a Schwarzschild black hole 



Dong-Hoon Kim 

Department of Physics, PO Box 118440, 
University of Florida, Gainesville, FL 32611-8440 
(Dated: February 7, 2008) 

Abstract 

The interaction of a charged particle with its own field results in the self-force on the particle, 
which includes but is more general than the radiation reaction force. In the vicinity of the particle 
in curved spacetime, one may follow Dirac and split the retarded field of the particle into two parts, 
(1) the singular source field, ~ q/r, and (2) the regular remainder field. The singular source field 
exerts no force on the particle, and the self-force is entirely caused by the regular remainder. We 
describe an elementary multipole decomposition of the singular source field which is an important 
step in the calculation of the self-force on a scalar-charged particle orbiting a Schwarzschild black 
hole. 

PACS numbers: 
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I. INTRODUCTION 



According to the equivalence principle in general relativity, a particle of infinitesimal mass 
orbits a black hole of large mass along a geodesic worldline F in the background spacetime 
determined by the large mass alone. For a particle of small but finite mass, the orbit is 
no longer a geodesic in the background of the large mass because the particle perturbs the 
spacetime geometry. This perturbation due to the presence of the smaller mass modifies 
the orbit of the particle from an original geodesic in the background. The difference of the 
actual orbit from a geodesic in the background is said to result from the interaction of the 
moving particle with its own gravitational field, which is called a self-force [1]. 

Historically, Dirac [2] first gave the analysis of the self-force for the electromagnetic field 
of a particle in flat spacetime. He was able to approach the problem in a perturbative 
scheme by allowing the particle's size to remain finite and invoking the conservation of the 
stress-energy tensor inside a narrow world tube surrounding the particle's worldline. Dcwitt 
and Brehme [3] extended Dirac's problem to curved spacetime. Mino, Sasaki, and Tanaka 
[4] generalized it for the gravitational field self-force. Quinn and Wald [5] and Quinn [6] 
worked out similar schemes for the gravitational, electromagnetic, and scalar field self-forces 
by taking an axiomatic approach. 

In Dirac's [2] flat spacetime problem, the retarded fleld is decomposed into two parts: (z) 
The flrst part is the "mean of the advanced and retarded flelds" which is a solution of the 
inhomogeneous fleld equation resembling the Coulomb q/r piece of the scalar potential near 
the particle, (ii) The second part is a "radiation" field which is a homogeneous solution of 
Maxwell's equations. Dirac describes the self-force as the interaction of the particle with 
the radiation field, a well-defined solution of the vacuum field equations. 

In the analyses of the self-force in curved spacetime, the Hadamard form of Green's 
function [3] is employed to describe the retarded field of the particle. Traditionally, taking 
the scalar field case for example, the retarded Green's function G"^^*(p,p') is divided into 
direct and tail parts: {i) The direct part has support only on the past null cone of the field 
point p. (ii) The tail part has support inside the past null cone due to the presence of the 
curvature of spacetime. Accordingly, the self-force on the particle consists of two pieces: (i) 
The first piece comes from the direct part of the field and the acceleration of the worldline 
in the background geometry; this corresponds to Abraham-Lorentz-Dirac (ALD) force in 
fiat spacetime. (ii) The second piece comes from the tail part of the field and is present 
in curved spacetime. Thus, the description of the self-force in curved spacetime reduces to 
Dirac's result in the fiat spacetime limit. In this approach, the self-force is considered to 
result via 



The quantities u, v and a are familiar from the Hadamard expansion of a Green's function 



Although this traditional approach provides adequate methods to compute the self-force, 
it does not share the physical simplicity of Dirac's analysis where the force is described 



from the interaction of the particle with the quantity [1] 

^self ^ ^ret _ ^dir 



(1) 





[3]. 
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entirely in terms of an identifiable, vacuum solution of the field equations: unlike Dirac's 
radiation field, the ip^'^^^ in Eq. (2) is not a solution of the vacuum field equation V^ip = 0. 
Moreover, the integral term in tjj^'^^^ comes from the tail part of the Green's function and 
is generally not differentiable on the worldline if the Ricci scalar of the background is not 
zero (similarly, the electromagnetic potential A^^'^ and the gravitational metric perturbation 
/i*^'' are not differentiable at the point of the particle unless (^Rab — \gabR^ and Rcadbu'^u'^, 
respectively are zero in the background [7]). Thus, some version of averaging process must 
be invoked to make sense of the self-force. 

In this paper we use an alternative method to split the retarded field in curved 
spacetime which is similar to Dirac's and suggested by Ref. [1]: (i) The singular source 
field ijj^ is an inhomogeneous field similar to the tidally distorted Coulomb field and exerts 
no force on the particle, (ii) The regular remainder field ip^ is a homogeneous solution of 
the field equation, analogous to Dirac's radiation field, and is entirely responsible for the 
self-force. This alternative split is reviewed briefiy in Section II. 

In Section III we give a brief overview of the mode-sum decomposition scheme to evaluate 
the self-force [8]. We consider a particle with a scalar charge q in general motion about 
a Schwarzschild black hole. A spherical harmonic decomposition provides the multipole 
components of of both ip^^ and tp^. Then, the mode by mode sum of the difference of these 
components determines and, thence, the self-force. The multipole components of ■0'^^* can 
be determined numerically while the multipole components of ip^ are derived analytically. 
In particular, the multipole moments of are generically referred to as regularization 
parameters [8]. This paper focuses on the analytical task of finding these regularization 
parameters. Our analytical results are summarized at the end of the Section in Eqs. (12)- 
(19). These results are in agreement with those of Barack, Ori, Mino, Nakano, and Sasaki 
[9, 10, 11]. 

The description of ip^ becomes particularly simple in a specially chosen co-moving frame: 
the THZ normal coordinates, introduced by Thorne and Hartle [12] and extended by Zhang 
[13], are locally inertial on a geodesic. In Section IV we obtain a simple form for tp^ us- 
ing the THZ coordinates and then re-express it in terms of the background Schwarzschild 

coordinates. 

Section V outlines our derivation of the regularization parameters which, while not ele- 
mentary, appears to us to be more compact than the derivations of others [9, 10]. 

Appendix A provides some mathematical details concerning the hypergeometric functions 
and the different representations of the regularization parameters in connection with them. 

Notation: {t, r, 9, (p) are the usual Schwarzschild coordinates. The particle moves along 
a worldline F, parameterized by the proper time r. The points p and p' refer to a field point 
and a source point, respectively, on the worldline F of the particle. In the coincidence 
limit p p'. The coordinates (T, X, Y, Z) are intermediate coordinates derived from the 
Schwarzschild coordinates, while (T, A", Z) are the THZ coordinates attached to the 
worldline V of the particle with p = V'^^ + + 2,'^- 

II. DECOMPOSITION OF THE RETARDED FIELD 

The recent analysis of the Green's function decomposition by Detweiler and Whiting [1] 
shows a method to split the retarded field into two parts 

^-•^t = V'^ + V'^, (3) 
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where and ip are the singular source field and the regular remainder field, respectively. 
The source function for a point particle on the worldline F is 

Q{p)^qj{-g)-'/'5\p-p'{T'))dT'. (4) 

The singular source field is an inhomogeneous solution of the scalar field equation 

= -47r^ (5) 

in the neighborhood of the particle, just as ijj'^'^^ is. And is determined in the neighborhood 
of the particle's worldline entirely by local analysis, t/'^i defined by Eq. (3), is then necessarily 
a homogeneous solution and is therefore expected to be differentiable on T. According to 
Ref . [1] , will formally give the correct self-force when substituted on the right hand side 
of Eq. (1) in place of t/)®^'^. In this paper we adopt this decomposition, and determine an 
analytical approximation, via a multipole expansion, of ■0^, which is to be subtracted from 
^ret explicit computation of the self-force. 



III. MODE-SUM DECOMPOSITION AND REGULARIZATION PARAMETERS 



By Eq. (1) the self-force can be formally evaluated from 



■self 



lim 

q hm Va('0''''* - V'^ 



(6) 



where p' is the event on F where the self-force is to be determined and p is an event in the 
neighborhood of p'. For use of this equation, both T'^^^ip) and ^aip) would be expanded 
into multipole £-modes, with J-'fa^{p) determined numerically. 

For the Schwarzschild spacetime, the source function g{p) is expanded in terms of spher- 
ical harmonics, and a similar expansion for is 



ret 



im 



(7) 



where t) is found numerically. The individual components in this expansion are 

finite at the location of the particle even though their sum is singular. Then, the £ component 



JF^^'f is finite 



em- 



(8) 



The singular source field is determined analytically in the neighborhood of the particle's 
worldline via local analysis (see Section IV) and its mode-sum decomposition provides 

^fa-Q'^aT.i^LYem, (9) 

m 

which is also finite at the location of the particle. Eqs. (6), (8) and (9) now imply that 

TT'^ = T.^^,[^iTip)-^lip) 

= q Y: lim V„ Y.i^Z - i^lJyem (10) 



p^p' 
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evaluated at the location of the particle. 

We follow Barack and Ori [8] in defining the regularization parameters, except that the 
singular source field ip^ is used in place of ip'^^^, 



(11) 



In Section V these regularization parameters are derived from the multipole components of 
evaluated at the source point, 



Ar = -sgn(A) 



r2l + J2/r2' 



"3/2 



1 + jyri 



3F. 



5/2 



Br 



Fl/2 



5/2 



2 (1 + jyrl) 



2^5/2 



«gn(A)!^ . , (12) 

(13) 

0, (14) 

(15) 

+ ^ 1— W + ^ , . .,./o' ) , (16) 

.(1 + J'lrlf^" ^ 2a + J7^. ' ^^^^ 

a = a = = 0, (18) 

Ag^Bg^Ce^ 0, (19) 

where /S. = r — E = —Ut = (1 — IMjr^ {dt/dr)^ (r: proper time) and J = — 
r1 {dcp/dr)^ are the conserved energy and angular momentum, respectively, and r = — 
(dr/dr)^. The subscript o denotes evaluation at the location of the particle. Also, shorthand 
for the hypergeomctric function is Fp = 2F1 (^p, |; 1; J^/(r^ + J^)^ (see Appendix A for more 
details aboiit the hypergeomctric functions and the representations of the regularization 
parameters in terms of them). 



2 (1 + jyr^.f^ 



Fi/2 — Fz/2 3(F5/2 — F3/2) 



IV. DETERMINATION OF IN LOCALLY INERTIAL COORDINATES 

In the vicinity of an event p' on a timelike worldline F, physics is most easily described 
in terms of locally inertial coordinates, where the time coordinate T on F is equal to the 
the proper time, and the orthogonal Cartesian-like spatial coordinates are (A*, 3^, Z) and 
centered on F. At p', with locally inertial coordinates the spacetime metric is just the flat 
Minkowski metric, and all of the Christoffel symbols vanish. Locally inertial coordinates 
are not unique and have an ambiguity at O(p^), where = \/ + y"^ + Z"^. For example, 
differences of 0{p^) distinguish Riemann normal from Fermi normal coordinates [14]. For 
our purposes a locally inertial coordinate system introduced by Thorne and Hartle [12] 
and later extended by Zhang [13] is particularly advantageous. In these THZ coordinates 
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Detweiler, Mcssaritaki, and Whiting [7] (cited henceforth as Paper I) show that the scalar 
wave equation takes a simple form and also that 



V'' = g/p + o(pV7^^), 



(20) 



where TZ represents a length scale of the background geometry (the smallest of the radius 
of curvature, the scale of inhomogeneities and time scale for changes in curvature along V). 
Approximation (20) is accurate enough for self-force regularization because 



v,^s^ v.(g/p) + o(p/7^3), 



(21) 



and the 0{p/7i?) remainder vanishes at the particle. 

For the derivation of the regularization parameters from the multipole components of 
V a'^^i requires that p in Eq. (20) be expressed in terms of the coordinates of the background 
geometry. Thus, we look for the relationship between the original Schwarzschild coordinates 
{t, r, 9, (f)) and the THZ coordinates (T, A", Z) associated with an event p' on F. 
However, in Section V we surprisingly find that any locally inertial coordinate system is 
sufficient to determine the regularization parameters quoted in Eqs. (12)-(19), which are 
actually independent of the 0{p^) terms in the THZ coordinates. An elementary discussion 
in Weinberg [15] determines this coordinate transformation through terms of 0{p^) in two 
steps: 



(i) Find inertial coordinates in the neighborhood of the event p' on F in terms 
of a Taylor expansion of the Schwarzschild coordinates about p', where the Schwarzschild 
coordinates at p' are and the subscript o denotes evaluation at p'. Weinberg's [15] 
Eq. (3.2.12) is 



1 



= x^ + M^„(x« - x^) + -M^a nx K - - + o[{x - x,n 



(22) 



where we may choose X^ — and M^^ = diag 

as this choice recenters and rescales the Schwarzschild coordinates to T 
X = M^,(r - To), Y = M^^{(t> - (f),), Z = M^eie - 9,). 



for convenience 



(ii) Boost with w^, the particle's four- velocity at p' as measured in this Carte- 
sian frame, to obtain the final coordinates ; 



K^'aX^ 



+ 0[{x-x,fl (23) 



where 



A^ 



u 



sym 



(24) 



is the upper half of the symmetric matrix A^' a with u'^ = (u^y + ('?^^)^ + {u^^ [16]. 
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With the choice of M^t = (1 - ^M/r^f^, M^r = (1 - 2M/ro)"^/^ M^^ = To sin and 



M^g = -ro, it follows that 



so that 



dg 



A'B' 



(25) 



(26) 



Eqs. (25) and (26) arc the desired locally inertial features for a particle in the Schwarzschild 
geometry at event x'^ with four- velocity u°'. 

To simphfy the calculations, we confine the particle's orbit to the equatorial plane 9o — 
7r/2 and have 



= diag /-V2^ 



(27) 



where / = (1 — ^). This constraint to the equatorial plane makes u — and we rewrite 

^ (28) 



^ {u^, u') = o) , 



in terms of the Schwarzschild coordinates and the constants of motion: E = —Ut ~ 
fidt/dr)^ and J = — rl{d(f)/dT)^ are the conserved energy and angular momentum, 
respectively, and f = — (dr/dr)^. Prom this it follows that — f~^E'^ — 1 and we have 



l + rV(/V2E + /) Jr/[r,{E + f/^)] 
l + jy[rlir'/'E + l)] 
sym 1 



(29) 



Now we are able to express in terms of the Schwarzschild coordinates using Eq. (23) 
and obtain 



^ X'Xj = 6ijA'cA\M'^,M''d [{x'' - xl){x^ - <) + r^,|„ {x^ - 0(a;^ - x'J{x'' - 

+0[{x-x,)% (30) 

where /, J — 1, 2, 3. Then, after a substitution from Eqs. (27) and (29), becomes 

2Er 

= {E^ -m- to)' -^{t- Q{r - To) - 2EJ{t - Q{4> - 4>o) 



f 



2Jr 



+ ^[^ + ^]ir-ror + ^{r-r,){(l>-<Po) + {K + J'){<P-Ay + K[e-- 



f\ f 



f 



n 



MEf , ,o M/ 2^2 r^x 



o 

MEr 



^r^Ef{t - tM - <^of + r^Ef{t - Q (^^ - - 



{t - Q{r - rM - (f),) 
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-(r - Tofi^ - (f>o) 

+ 0[(x-Xo)"]. 



TT 



2 



) 



2 




(31) 



The substitution of Eq. (31) into Eq. (20) approximates in terms of the Schwarzschild 
coordinates and leads to the derivation of the regularization parameters in the next section. 

In the above analysis the 0[{x — Xo)^] term in Eq. (23) contributes to the 0[{x — Xo^] 
terms of Eqs. (30) and (31). To the level of accuracy we desire for the mode-sum 

regularization parameters in this paper, that is to say, to the determination of C^-terms, it 
is not necessary to specify the 0[{x — Xo)^] term in and, hence, not necessary to specify 
the 0[{x — Xo)^] terms in the spatial THZ coordinates X, y, Z. 

V. REGULARIZATION PARAMETERS FOR A GENERAL ORBIT OF THE 
SCHWARZSCHILD GEOMETRY 

In Section IV, we have seen that an approximation to ip^ is 



Following Paper I [7], the regularization parameters can be determined from evaluating the 
multipole components of da{q/p) {a = t, r, 9, (j) for the Schwarzschild background). The 
error, O^p^ /IZ^) in the above approximation is disregarded since it gives no contribution to 
VaV'^ as we take the "coincidence limit" , — > , where x"- denotes a point in the vicinity 
of the particle and a;" the location of the particle in the Schwarzschild geometry. 

In evaluating the multipole components of da{q/p), singularities are expected with certain 
terms. To help identify those singularities, we introduce an order parameter e which is to 
be set to unity at the end of the calculation: we attach e" to each 0[{x — Xq)"] part of in 
Eq. (31) and re-express as 



where Vu, Vm, and Vw represent the quadratic, cubic and quartic order parts of p^, re- 
spectively. Here we pretend that the quartic part Viy is also specified: this will help us 
to perform the structure analysis for Ca-terms later in Subsection VC when we prove that 
these regularization parameters always vanish. 

We express da{l/p) in a Laurent series expansion where every denominator of this expan- 
sion takes the form of Vi/"^ (n = 3, 5, 7). Thus, Vn plays an important role in the multipole 
decomposition, and the quadratic part Vu, directly taken from Eq. (31), is not yet fully 
ready for this task. First, — 0o must be decoupled from r — Tq so that each appears only 
as an independent complete square. Coupling between t — and — 0o docs not create 
difficulty in the decomposition. Thus, we reshape the quadratic term of Eq. (31) into 



V'' = g/p + 0(pV7e3). 



(32) 



(33) 



2Errl 



{t-QA-2EJ{t-Q{(l>-cf>') 



f {r^o + J') 




TT 



2 



) 



2 



(34) 



8 



with 

^ '^o - J^^y (35) 

where A = r — Tq, and an identity = — f {1 + J^/rl) is used for simphfying the 
coefficient of A^. Here, taking the coincidence hmit A — > 0, we have 0' — > 0o- This same 
idea is found in Mino, Nakano, and Sasaki [9]. Also, for the multipole decomposition the 
quadratic part must be analytic and smooth over the entire two-sphere, and we write 

Vii = {E^ -m- Uf - . , 2 ^ °72N - ^o)A - 2EJ{t - to) singsin(0 - ct>') 



+0[{x-x,f]. (36) 

Here we have used the elementary approximations (j) — cp' = sin(0 — cp') + O[(0 — cp'Y] and 
1 = sin^ + 0[(^-7r/2)2]. 

To aid in the multipole decomposition we rotate the usual Schwarzschild coordinates 
by following the approach of Barack and Ori [10] such that the coordinate location of the 
particle is moved from the equatorial plane ^ = 7r/2 to the new polar axis. The new angles 
© and $ defined in terms of the usual Schwarzschild angles are 

sin 6 cos(0 — 0') = cos O 
sin^sin(0 — 0') = sin© cos $ 

cos 6' = sin© sin $. (37) 

Also, under this coordinate rotation, a spherical harmonic Y(^^{0, (p) becomes 

YUO,<P)= E ai^,Y,m'ie,<^), (38) 

m'=-e 

where the coefficients a^^/ depend on the rotation {9, (p) (©, $) as well as on m and 
m', and the index ^ is preserved under the rotation [17]. As recognized in Ref. [10], there is 
a great advantage of using the rotated angles (©, $): after expanding da{q/ p) into a sum of 
spherical harmonic components, we take the coincidence limit A — > 0, © ^ 0. Then, finally 
only the m = components contribute to the self-force at © = because Yem{0, $) = 
for m 7^ 0. Thus, the regularization parameters of Eq. (11) are just {£, m — 0) spherical 
harmonic components of da{q/p) evaluated at x'^. 

Now, with these rotated angles, Vn is re-expressed as 

Vu = {E'-f){t-toY- r, , °M -QA - 2E Jit- to)smecos^ 

J [r^ + J ) 



^ / 9 sin^ $ ' 



2/2 (^2 + J2)2 



r^E^A^ 

+ 1-C0S© 



J2 sin^ <J> 



+0[ix-x,)% (39) 
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where the elementary approximation sin^ © = 2(1 — cos©) + 0(©^) is used. We may now 



define 



p2 = ^E'_f^^t-Q'- 
+2 {rl + J') (l 



° (t-to)A-2EJ(t-to)sin©cos$ 



J2 sin^<l>^ 
r2 + J2 J 



2/^ (r2 + J2)^ (l - 



+ 1 - cos © 



In particular, when fixing t = to, we define 



t=to 



with 



= 2(r^ + j2)x((^' + l-cos©) 
J2 sin^ $ 



rl + J-^ 



and 



5^ = 



2/2 (r2 + J2)2^- 

Now we rewrite Eq. (33) by replacing the original quadratic part Vu with p 



;;2 



(40) 

(41) 
(42) 
(43) 

(44) 



where Viy now includes the additional quartic order terms that have resulted from the 
replacement of Vu by p^. A Laurent series expansion of da{'\-/ p)\^^^^ is 



t=to 



Pi 



\t=to ^-2 I 



Pi 



Pi 



^|e-^ + 0(e°). (45) 



After the derivatives in Eq. (45) are taken, the dependence upon ©, $ and r may be 
removed in favor of po, X ^ by use of Eqs. (41)-(43). Then the three steps of (i) a 
Legendre polynomial expansion for the dependence, while r and $ are held fixed, followed 
by {ii) an integration over while r is held fixed, and finally [iii) taking the limit 5 — > 0, 
together provide the regularization parameters. The techniques involved in the Legendre 
polynomial expansions and the integration over $ are described in detail in Appendices C 
and D of Paper 1 [7] . 

Below in Subsections VA and VB, we present the key steps of calculating the Aa and 
Ba regularization parameters in Eqs. (12)-(19). 



A. ^a-terms 

We take the term from Eq. (45) and define 

^ Po 

Then, we proceed with our calculations of the regularization parameters one component at 
a time. 
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1. At-term: 

First we complete the expression for Qt[e~^] by recalling Eqs. (40) and (41) 



-\Po' {~p 



t=to 



2 



q^ErrlAx (^2^. r>\"^/^ 

( d) + 1 — cosGj 



-3/2 / 2Errl/^ 



+ 2£;j sine cos $ 



2V2/ (r2 + J2)5/2 

q^EJx'^^'^ cos $ (9 
~ V2(r2 + J2)3/2 M 



(5^ + 1 - cose) 



-1/2 



(47) 



where c)/(9e|^ means that A is held constant while the differentiation is performed with 
respect to e. 

According to Appendix D of Paper I [7] , for p > 1 



(5^ + 1 - cose) 
and for p = 



-P-i/2 ^ ~ 2£+ 1 



(52 + l-cose) = ^[\/2 + C>(£5)]P^(cose), 5^0 



[1 + 0{U)\ P^(cose), 5^0, (48) 



^=0 



(49) 



Then, by Eqs. (48) forp = 1, (49) and (43), in the limit 5^0 (equivalently A ^ 0) Eq. (47) 
becomes 



limg,[e-2] = sgn(A) 

(r2 + J2) / 



{rl + J2) 



3/2 



g(£+-jP,(cose) 

p^cose). 



(50) 



Then, we integrate limA^oQt[e ^] over $ and divide it by 27r (we denote this process by 
the angle brackets "()") 



(lim = sgn(A)i|r^| + 1) P,(cose), 



(61) 



where we exploit the fact that ^^"^ cos = to get rid of the second term in Eq. (50) 
[19]. Appendix C of Paper I [7] provides (x"^) = 2F1 (l, \] 1; a) = Pi = (1 - Q;)~^/^ where 
a = J^/ {rl + J^). Substituting this into Eq. (51), the regularization parameter At is the 
coefficient of the sum on the right hand side in the coincidence limit e — > 



At = sgn(A) 



r2l + J7r2- 



(52) 
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2. Ar-term: 



Similarly, we have 



t=to 



(53) 



Here, before computing dr (p^)!^^^^ we reverse the steps of Eqs. (34), (36), (39) and (40) to 
obtain the relation 

p' = Vn + 0[{x-x,)% (54) 

where Vu is now back to Eq. (34). Differentiating this with respect to r and going through 
the steps of Eqs. (36) and (37), Eq. (53) can be expressed with the help of Eq. (41) as 



Qr[€~ 



p 



1 — cos 9 



-3/2 



— + fJr sin 9 cos $ 

ri + 



(55) 



Then, the rest of the calculation is carried out in the same fashion as for the case of At-texvn 
above. We obtain 



A = -sgn(A) 



E 



rlf{l + J^/riy 



(56) 



3. Aj,-term: 



First we have 



t=to 



(57) 



Taking the same steps as used for A^-term above via Eqs. (54), (36) and (37) in order, we 
obtain 

(f) = 2 (r^ + J2) sin 9 cos $ + 0[{x - x^f] . (58) 

Then, in a similar manner to that employed in the previous cases, in the limit A — > 
Eq. (57) becomes 



P^(cos9). 



(59) 



The right hand side vanishes through "()" process because {x ^'^^cos$y = 0. Hence, 

= 0. (60) 



4- Ag-term: 

It is evident from the particle's motion, which is confined to the equatorial plane 9o — ^, 
that no self-force is acting on the particle in the direction perpendicular to this plane. This 
is due to the fact that both the derivatives of retarded field and the singular source field 
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with respect to 9 tend to zero in the coincidence hmit. Our calculation of Aq should support 
this. Through the same process as employed before, we have 



t=to 



with 



t=to 



2rl sinQsin^ + 0[{x - Xof]. 



P^(cose). 



(61) 
(62) 

(63) 



Then, similarly as in the case of A^-ierva. above 

Again, via "( )" process, the right hand side vanishes because {xT^^'^ sin $^ = 0. Thus, 

Ae = 0. (64) 



B. Ba-terms 



We take the e ^ term from Eq. (45) and define 

where for computing da (p^), Eq. (54) should be referred to, and Vm is the cubic part taken 
directly from Eq. (31). 

We may express this in a generic form 

Qa[e-'] = E E E (66) 

n=l fc=0 p=0 Po 

where A = r — Tq, and bnkp(a) is the coefficient of each individual term that depends on n, k 
and p as well as a, with a dimension TZ''~^ for a — t,r and TZ^ for a — 9, cf). We recall from 
Eqs. (34) and (35) that the first of the steps to lead to in the denominator is replacing 
— 0o by (0 — (f)') — JrA/ f{rl + J"^) to eliminate the coupling term A (0 — 0o)- This makes a 
sum of independent square forms of each of A and — 4>', which is a necessary step to induce 
the Legendre polynomial expansions later. Thus, to be consistent with this modification in 
the denominator, — 0o in the numerator on the right hand side of Eq. (66) should be also 
replaced by (0 — 0') — JrA/f{rl + J^). Then, this will create a number of additional terms 
apart from (0 — 0')™" when we expand the quantity [(0 — 0') — JrA/f{r'^ + J^)] raised, say, 
to the m-th power, and the computation will be very complicated. 

By analyzing the structure of the quantity on the right hand side of Eq. (66) one can 
prove that — 0o may be replaced just by — 0' in the numerator without the term 
—JrA/f{r1 + ,P) (the same idea is found in Mino, Nakano, and Sasaki [9]). The verification 
follows. The behavior of the quantity on the right hand side of Eq. (66), according to the 
powers of each factor, is 

Qa[e-'] ~ p-(2-+i)A2'^-^ (0 - 4>,f-'' {9 - 0''7^^ (67) 
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where s — k — 1 ior a — t, r and s — k ior a — 9, (f). Further, 

-| k-2p 

(0 - ^Po)'-'" 



(0-00 ''^ 



k-2p 

= E Cfcpi (0 - 00' A'^-^f- ~ (0 - 0')^ A^-2f-V7^'=-2^'-' (68) 



i=0 



(sin 0)' (cos /^k-2p-i i^k-2p-i j^o[{x- Xo)'^~'^+'], (69) 



where a binomial expansion over the index i = 0, 1, ■ ■ ■ , A; — 2j9 is assumed with 
^^^fc_2p-i ^gg^)^ (^gg^) _ ^ry replaced by [sin(0 - 0')]' + 0[{(j) - 0')*^^] 

— the term 0[{x—Xo)'^~^^'^^] at the end results from this O[(0 — 00*^^] , then the coordinates 
are rotated using the definition of new angles by Eq. (37). Also, by Eq. (37) again 

(^e -^y' = (sin 6)'^ (sin*)'^ + 0[{x - Xo)'^+']. (70) 
Using Eqs. (69) and (70), the behavior of Q[e~^] in Eq. (67) looks like 

Qa[e-'] - p-(^n+l)^2n-2p-^ ^^^^Q^P+i (^OS $)' (sin$)'^7^^ (71) 

where s = 2p + i — 1 for a = t, r and s = 2p + i for a = 6, (j), and any contributions from 
0[{x — Xo)^~'^^^'^] in Eq. (69) and from 0[{x — sJq)^^^^] in Eq. (70) have been disregarded: 
by putting these pieces into Eq. (67) we simply obtain e^-terms, which would correspond to 
0{£~'^) in Eq. (11) and should vanish when summed over £ in our final self-force calculation 
by Eq. (10). Qafe"^] then can be categorized into the following cases: 

(i) i = 2j + l(j = 0,l,2,.--) 

The integrand for "()" process, F($) = (cos^)^^"*"^ (sin$)^^ has the property F($ + 
tt) = -F($). Thus 

[Qa[e-'])-0, (72) 



(ii) i = 2i(j = 0, 1, 2, •••) 

Using Eqs. (41) and (43), we can express (sin©)^*'"''* in Eq. (71) above in terms of po 
and A via a binomial expansion 

(sin efP+^' = [2(1- cos 0)]^+^' + 0[ix - Xo)'(^+^')+'] 

= E d,,,pl'^A'^P+^-'^^ + 0[{x - Xo)2(^+^''+2] (73) 

q=0 

where g = 0, 1, ■ ■ ■ , p + j is the index for a binomial expansion and dpjq ~ 1/TZ'^^p~^^\ 
When Eq. (74) is substituted into Eq. (71), the contribution from 0[{x - Xo)^^+^')+^] 
can be disregarded since it would correspond to O(e^) again. Then, we have 

QJe-^] ~ (sin <l>)'^ (cos $)'^' p-2{«-g)-i/^2(n-g)^.^ ^^^^ 

where s = — 1 for a = t, r and s = for a = 6, (f), and we can guarantee that n — q>0 
always since 0<q<p + j= p+^i, 0<i<k — 2p and p < k < 2n. Then, Eq. (75) 
can be subcategorized into the following two cases; 
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(ii-1) n — q > 1 

By Eqs. (41), (43) and (48) 



Qa[e~^] ~ (sin (cos $)^^ AP^(cos e)7^" — > 0, (76) 



(ii-2) n-q = 

By Eqs. (41), (43) and (49) 



Qa[e"'] ~ (sin$)'''(cos$) ■'P^(cose)7^^ (77) 

where s = — 1 for a = t, r and s = for a = ^, 0. 

Therefore, by analyzing the structure of (5a[c~^] we find that the e~^-tcrms vanish in all the 
cases except when n — q = 0. The non- vanishing Ba-terms are derived only from this case. 
Then, hyO<q<p + j= p + ^i,0<i<k — 2p and p < k < 2n together with n = q one 
can show that 

< k-2p-iandk-2p-i<0, i.e. k - 2p - i ^ 0. (78) 

Substituting this result into Eq. (68), then into Eq. (66) we may conclude that in the 
numerator of (5[e~^] in Eq. (66) one can simply substitute 

(0 - 0o)'"'" ^ (0 - 0')'"'" • Q- E. D. (79) 

The significance of this proof does not lie in the result given by Eq. (79) only, but also in 
the fact that the non- vanishing contribution comes only from the case n — q for Eq. (75), 

Qa[e-'] ~ (sin$)2^ (cos$)'("-^^ p-'n% (80) 

where n — 1, 2 and < p < n, and s = — 1 for a = r and s = for a = ^, 0. 

Below are presented the calculations of Sa-terms of the regularization parameters by 
component, in a similar manner to those for A^-terms. 



1. Bt-term: 
We begin with 

(81) 



t=to 



The subsequent computation will be very lengthy and it will be reasonable to split Qt[e ^] 
into two parts. First, let 

Qt(i)[e-']^-^~p-'dtVui\t=t., (82) 
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where 

dtVm\ 



t=to 



-2 1 



M\ EJA 



To J fro 



{4> - 4>o) + roEr 



(83) 



As proved at the beginning of this Subsection, every (0 — ^o)™ in the numerators of the e~^- 
term can be replaced by (0 — 0')"* without affecting the rest of calculation. Then, followed 
by the rotation of the coordinates via Eq. (37) 



MErA^ 






(-9 



fro 



sin © cos $ + 2roEf (1 - cos ©) 



+0 



{X - Xq) 

~pl 



(84) 



where an approximation sin^ = 2(1 — cos G) + 0[{x — XoY] is used to obtain the last term 
inside the first bracket. Here we may drop off the term 0[{x — XoY / Po] , which is essentially 
O(e^), for the same reason as explained at the beginning of this subsection. Then, using the 
same techniques as used to find Aa-terms, we can reduce Eq. (84) to 



q^MEf q^rlE^rx ^ 



2/2 {rl + J^r 



v^/ro {rl + J'f 
q^ErroX~^ i 



2(ro' + J')x(<^' + l-cose)]"''" 

(52 + 1-cose) 



2(r2 + J2)' 



(85) 



As we have seen before, by Eq. (48) (5^ + 1 — cos 0) ~ A ^ in the limit A ^ and the 
first term on the right hand side will vanish. The second term will also give no contribution 
to the regularization parameters because {xT"^^"^ cos $^ = 0. Only the last term, which is 
~ will give non-zero contribution according to the argument in the analysis presented 
above (see Eq. (80)). Using Eq. (49) in the limit A ^ and taking "()" process, Eq. (85) 
becomes 



1 ^2 Erlx'^^"^) °° 
4^ LT^Peicose) 



(86) 



The identity (x"*') = ( (l - ^sin^ ^) = 2F1 (p, |; 1, a) = Fp, with a = J^/ {rl + J^) is 



taken from Appendix C of Paper I [7], and we take the limit — > 



ErF- 



3/2 



2r2(l + j2/^2)3/2- 



(87) 



Now the remaining part is 



Qti2)[e-']^^Po' 



III 



t=to 



(88) 



where 
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III 



t=to 



2EfA 



f 



2EJ (0 - 0o) 



X 



1 + 



MA3 



+ 2 



5M\ JrA2 



(0-0o) 



+ 0[(a;-a;o)% (89) 



Taking similar procedures as above, the non- vanishing contributions turn out to be 



lini(5t(2)[e = ^hni -g^£;j^rroPo^cos^$sin^© 



\A-.o2rol + JVr2 1 



X 



-2 



Ef 



2r2(i + j2/^2 



2^3/2 



-3/2 



5/2 



X 



1 + 77^5/ g/'*^"^^'- <^'" 



where all other terms than ~ Po ^ again have been dropped off during the procedure since 
they vanish either in the limit A ^ or through the "()" process. Then, using the identity 
{X'f) = 2F1 (p, i; 1, a) = Fp, we have 



(limg,(2)[6-^] 



3g2 



Ef 



'5/2 



e^O 2r„2(l + j2/^2)3/2 1^ 3/2 l + j2/^2^ 

By combining Eqs. (87) and (91), we finally obtain 



-Ef 



'3/2 



SE 



5/2 



2N5/2 



(91) 



(92) 



2. Br -term: 

Prom Eq. (65) we start with 

Qr[e-'] = q' 



2 53 4 -5 



PS 4 J 

Then, following the same steps as taken for the case of E^-term above, we obtain 



'1/2 



{l + J^rlf^ 2(l + JVr2) 



+ 



(1 - 2/-V2)F3/2 , 3f-'PF,/2 



2^3/2 



+ 



2 (1 + jyr^o) 



2^5/2 



(93) 



(94) 
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3. B^-term: 



Again, from Eq. (65) 



t=to 



(95) 



Then, similarly we can derive 



Bj, = — r 



Fi/2 — -F3/2 ^ 3(^5/2 — -F3/2) 



(96) 



Bg-term: 

As 745/ vanishes, so should Bg. Prom 



one finds that there is no term like ~ p~^: all terms are either like ~ A^^/Pq""*"^ or like 
~ /\2n-i gjj^ Q cos$/Pq"+^ (n = 1, 2), which vanish in the limit A ^ or through the "( )" 
process. Thus 

Bg = 0. (98) 



C. Ca-terms 

We have mentioned before that Ca-tcrms, which originate from e'^-tcrm in Eq. (45), 
always vanish. This can be proved by analyzing the structure of e°-tcrm. First we specify 
the e^-order term for f?a(l/p)lt=to ^ Laurent series expansion and define 

n f 01 ^ 2 f l daVMt=t^ 3 (d^Vui) Vui\^=^^ + [da (p')] Piv|,=,„ 



2 pl A PI 

'lll\t=to \ 



15 [damvf 



16 PI 
Generically, this can be written as 



(99) 



3 2n+i [m c„,,(„)A2-+i-'= (0 - <p,f-'^ {e - f ) 



2p 



= E E E — ^stF — (100) 

n=l A;=0 p=0 A'o 

where A = r — Tq, and Cnkp{a) is the coefficient of each individual term that depends on n, 
k and p as well as a, with a dimension T^'^^^ for a = t, r and 7^'^^^ for a = 6, (f). 

The behavior of (3a[e°], according to the powers of each factor on the right hand side of 
Eq. (100), is 

Qa[e'] ~ p-(2"+i)A2"+i-'= (0 - cp^)'-'^ (9 - ^y'n% (101) 
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where s — k — 2 for a = t, r and s = k — 1 for a — 9, (f). Following the same procedure as in 
the beginning of Subsection VB, Eq. (101) becomes 

Qa[e°] ~ A2"+i-2p-i (sin efP+' (sin $)'^ (cos 7^^ (102) 

where a binomial expansion over the index i — 0, 1, ■ ■ ■ , k — 2p is assumed, and s — 2p+i — 2 
for a — t, r and s — 2p + i — 1 for a — 9, (f). Here we have disregarded any by-products 

like 0[{x - Xo)''-'^P+^] and 0[{x - XofP+^l, which originate from (0 - (f)o)''~'^^ and (9 - f )^^, 
respectively when we rotate the angles: by putting them back into Eq. (101) we simply 
obtain e^-terms, which would correspond to 0{£^'^) in Eq. (11) and should vanish when 
summed over i in our final self-force calculation by Eq. (10). Then, the rest of the argument 
is developed in the same way as in the beginning of Subsection VB: 

(i) ? = 2j + l(j = 0, 1, 2, ■■■) 

The integrand for "()" process, i^($) = (cos $)^"''''^ (sin$)^^ has the property F($ -|- 
tt) = -F($). Thus 

{Qa[e'])^0, (103) 

(ii) ^ = 2j(j = 0, 1, 2, •••) 

We have 

Q„[e°] ~ (sin $)'^ (cos <5f p-2(n-,)-1^2(n-g)+l^.^ (^04) 

where g = 0, 1, • • • , p -|- j is the index for a binomial expansion and s — —2 for 
a — t, r and s = — 1 for a = ^, 0. Here we can guarantee that n — q > — |, i.e. 

n — g = 0, 1, 2, ■ ■ ■ since 0<q<p + j= p + ^i,0<i<k — 2p and p < k < 2n + 1. 
Then, Eq. (104) can be sub categorized into the following two cases; 

(ii-1) n-q>l 

By Eqs. (41), (43) and (48) 

Qa[e^] ~ (sin $)^^ (cos A2p^(cos 0)71' — > 0, (105) 

(ii-2) n-q^O 

By Eqs. (41), (43) and (49) 

Q„[e°] ~ (sin (cos AP^(cos e)7^" — > 0, (106) 

where s = — 2 for a = r and s = — 1 for a = 9, (p. 

Clearly, in any cases the quantity (5a[e°] does not survive, therefore we can conclude that 
Ca-terms are always zero. Q. E. D. 

Also, this justifies the argument that we need not clarify the term 0[{x — Xo)^] in Eq. (23) 
and its contribution to p^, which is 0[{x — Xo)^] in Eqs. (30) and (31) in Section IV or 
Viv in Eqs. (33) and (44) in Section V: by the analysis of the generic structure given 
above, — | c^a^ivlt=t„ /Po or | [c^a (p^)] ^ivlt=t„ /Po would simply vanish in the coincidence 
limit X — > Xo, regardless of what Vw is. 
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VI. DISCUSSION 



Self-force analysis in curved spacetime relies upon the ability to divide the field of a point 
charge into two parts. One part is singular and exerts no net force on the charge itself. The 
remainder is a regular, smooth vacuum field and is entirely the cause of any self-force. We 
see, in this manuscript, that the singular field is adequately approximated by its Coulomb 
field in a coordinate system which is locally inertial and centered upon the charge. With 
this elementary approximation of the singular field, it is guaranteed that the remainder 
is at least differentiable and provides the correct self-force. For a charge moving in the 
Schwarzschild geometry, the multipole moments of the singular field are the regularization 
parameters which are necessary for computing the self-force from a multipole expansion, and 
these regularization parameters are calculated with a relatively elementary procedure. Our 
analysis agrees with that of others [9, 10] and appears to us to be the most straightforward 
calculation of these important parameters. 

Future work will use a higher order approximation for the singular field which will result 
in a more accurate and more differentiable approximation for the regular remainder. In 
practice, a higher order approximation of the self-force significantly speeds up convergence 
in a mode sum, as demonstrated in Paper I [7] . 

The simplicity of our methods should also make them useful for self-force analysis in the 
context of the Kerr geometry. 

Another avenue for future work involves calculating gravitational self-force regularization 
parameters. Thus far, pubhshed parameters [9, 10] focus upon the self-force, rather than 
upon the the metric perturbations themselves. However, it is clear that a gravitational self- 
force calculation of a gauge invariant quantity requires the metric perturbations explicitly 
as well. 

Future work will find the THZ coordinates to higher orders, and calculate higher order 
regularization parameters which will provide faster convergence of the £ sums and corre- 
spondingly more accurate results. 
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APPENDIX A: HYPERGEOMETRIC FUNCTIONS AND REPRESENTATIONS 
OF REGULARIZATION PARAMETERS 

In Section V we define 



X= 1 — Oi sin^ $ 



(Al) 



with 



J2 



(A2) 



a = 



20 



And we use 



X 



2 ri"^ 



n Jo 



1 - a sin^ 0?$ 



= 2F1 [p,-;l,a] = Fp. 



In particular, for the cases P — k p — —h we have the following representations 



and 



2 

11 \ 2 - 
'1/2 = 2Fi(-,-,l]aj = -K{a) 



F-i/2 = 2Fi(-^.^.l;«) =^^(«), 



(A3) 



(A4) 



(A5) 



where K{a) and E{a) are called complete elliptic integrals of the first and second kinds, 
respectively. 

If we take the derivative of F1/2 with respect to k = ^/a. via Eq. (A3), we obtain 



/2 



or using Eq. (A4) 

However, Ref. [18] shows that 



dk 
dk 



k ' 



dk 



OK 



E 



K 



dk k{l-k^) k 
Thus, by comparing Eq. (A7) and Eq. (A8) we find the representation 

2 E 2 E 



(A6) 



(A7) 



(AS) 



F3/2 — 



ttI — A;^ ttI — a 



(A9) 



Further, we can also find the representation for F^/2. First, taking the derivative of F3/2 
with respect to A; = ^/a via Eq. (A3) gives 



dF, 



3/2 



dk 



3-^3/2 _|_ 3F5/2 



k 



k 



(AlO) 



Also, using Eq. (A9) together with Eqs. (A3)-(A5), another expression for the same deriva- 
tive is obtained solely in terms of complete elliptic integrals 



dE 



3/2 



2{i + e)E -{i-k^)k 



dk IT k{l- A;2) 
Then, by Eqs. (A9), (AlO), and, (All) we find 



(All) 



F5/2 = — 



"2(2-a)^ k 



:i-a) 



1 — a 



(A12) 
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Now, using Eqs. (A4), (A9) and (A12), we may rewrite the non-zero Ba regularization 
parameters, Eqs. (15)-(17) in Section III as 



■? Ef 



K{a) - 2E{a] 



B, 



7r(l-2M/ro)(l + JVr2f/' 
ry2 f\k(a) - E(a] 



ro7r(J/ro)(l + JVr2)i/2' 
which are exactly the same to the results of Barack and Ori [10] . 
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